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Abstract 

\ We consider a version of the symmetric Anderson impurity model (compact- 



o 



ified) which has a non-Fermi hquid weak couphng regime. We find that in the 



\ Majorana fermion representation the perturbation theory can be conveniently 

' developed in terms of Pfaffian determinants and we use this formalism to cal- 

' culate the impurity free energy, self energies, and vertex functions. We also 

' derive expressions for the impurity and the local conduction electron charge 

and spin dynamical susceptibilities in terms of the impurity self-energies and 



' vertex functions. In the second-order perturbation theory, a linear tempera- 

Sin ■ 



ture dependence of the electrical resistivity is obtained, and the leading cor- 
rections to the impurity specific heat are found to behave as TlnT. The 

' impurity static susceptibilities have terms in InT to zero, first, and second or- 

! 

der, and corrections of In^ T to second order as well. The conduction electron 
static susceptibilities, and the singlet superconducting paired static suscepti- 
bility at the impurity site, are found to have second-order corrections InT, 
which we interpret as an indication that a singlet conduction electron pairing 
resonance forms at the Fermi level (the chemical potential). When the per- 
turbation theory is extended to third order logarithmic divergences are found 
in the only vertex function ro,i,2,3(0, 0, 0, 0), which is non-vanishing in the 
zero frequency limit. We use the multiplicative renormalization-group (RG) 
method to sum all the leading order logarithmic contributions. These give rise 



to a new weak-coupling low-temperature energy scale Tc = Aexp — 5 

which is the combination of the two independent coupling parameters. The 

RG scaling equation is also derived and shows that the dimensionless coupling 

constant U = ^ is increased as the high-energy scale A is reduced, so our 

perturbational results can be justified in the regime T > Tc. Below Tc the 

perturbation theory breaks down. 
75.15.Qm,71.45.-d,75.15.Nj 
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I. INTRODUCTION 



Strongly correlated electron systems, and especially the high cuprate superconductors, 
display features which have been difficult to reconcile with conventional theories of Fermi 
liquids and superconductivity. The various ideas that have been put forward to explain 
the behavior of these systems remain controversial. The unusual normal state of the high 
Tc materials has been interpreted as some form of non- Fermi liquid (non-FL), but there 
have been no fully convincing microscopic derivations of such a state in two dimensional 
systems with short range interactions. Understanding the normal state is rather important 
as it may be the key to understanding the nature and origin of the superconductivity. 
Anomalous behavior that has been observed experimentally in certain heavy fermion U- 
based superconductors has also been interpreted as some form of non-FL behavior One 
of the most striking characteristics of the non-FL in both these types of systems is the 
linear temperature dependence of the electrical resistivity. A phenomenological marginal- 
FL spectrum for the spin and charge fluctuations was put forward by Varma et al. [0 as a 
unified way of interpreting the diverse anomalies observed in the cuprate superconductors. 
The essential point in this theory is that the frequency dependence of the polarization is 
singular, while the momentum dependence is taken to be smooth. This picture is very 
similar to the two dimensional Luttinger liquid hypothesis used by Anderson |^. There is 
no generally accepted microscopic theory to explain the linear temperature dependence of 
the electrical resistivity. 

The lack of progress in developing a microscopic theory is due to the difficulty in solving 
the lattice models in strong coupling, and it may be quite a time before suitable methods 
can be developed to overcome these difficulties. There are impurity models, however, which 
display non-FL behavior; these are more accessible, and may provide valuable insights. 
For many of these models we have exact solutions Q], and in this context, the non-FL 
thermodynamic behavior of the two-channel Kondo model has been extensively investigated 
p|^IO|. However the exact solution for this model obtained from conformal field theory ^ 
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gives a resistivity that behaves as p{T) = p(0)(l + aVT) in the low temperature hmit, so 
it is not clear that the experimental observations can be interpreted using this theoretical 
model. The question then arises as to whether the linear temperature dependence of the 
resistivity, and the marginal-FL behavior, can be found in any single-impurity model. 

In this paper we examine a model, originally put forward by Coleman and Schofield ||10|| , 
which we will show displays just this type of behavior, at least in weak coupling perturbation 
theory. We can obtain the model starting from usual symmetric Anderson single-impurity 
model in the Majorana representation. The usual model has 0(4) symmetry arising from 
a product of the SU(2) spin and the SU(2) particle-hole symmetry. The modified model is 
obtained by reducing the 0(4) symmetry of the hybridization term to 0(3) symmetry. We 
then construct a systematic non-FL perturbation theory around the weak coupling limit. 
In contrast to the two-channel Kondo model, the linear temperature dependence of the 
electrical resistivity is obtained from second order perturbation theory. The result would 
appear to be valid at temperatures above a characteristic temperature Tc. At Tc, the sum 
of the leading order singular perturbational contributions to the electrical resistivity diverge 
indicating a breakdown of the perturbation theory at this energy scale. A brief report of 
this work was given in our previous paper |ri[| . 

The paper is arranged as follows. In section II, the modified Anderson single-impurity 
model is introduced and the general features of this model are discussed, in particular, there 
is an important relation between the conduction electron and impurity Green functions. In 
section III, we fully investigate the unperturbed Hamiltonian, and show that it displays 
a non-FL behavior in the spin and charge density- density spectrum. In section IV, we 
present the a new perturbation formalism in which the impurity free energy, single-particle 
Green functions, and two-particle vertex functions can be expressed in terms of Pfaffian 
determinants, a specific determinant defined from the ordinary antisymmetric determinant. 
In section V, using this formalism, expressions are derived for both impurity and local 
conduction electron spin and charge dynamical susceptibilities, including the conduction 
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electron singlet superconducting paired susceptibility, in terms of the impurity self energy 
and vertex functions. In section VI, the lower order results are given and discussed in detail. 
In section VII, multiplicative renormalization-group method is applied to this model to sum 
all the leading logarithmic contributions and obtain the scaling equation. Finally, some 
discussions and conclusions are given in section VIII. 



II. THE MODEL HAMILTONIAN 

The ordinary symmetric Anderson single-impurity model can be expressed in the form: 

H = itY, [Cl{n + l)a(n) - H.c] +iVY. [cl(0)rf. - H.c 

+U [d\d, - [d\d, - , (1) 

where the symmetric condition = —U /2 has been used, and the chemical potential is set 
to zero, the Fermi level. The Hamiltonian has 0(4) symmetry due to the SU{2) symmetry 
from the spin rotational invariance and an additional SU{2) from particle-hole symmetry, 
giving 0(4) ~ SU{2) ® SU{2). 

It is important to observe that the 0(4) symmetry can be explicitly displayed when the 



fermions of each type of spin are expressed in terms of four Majorana (real) fermions [jTO|: 



d^ = -^{di-id2), di = -^{-ds-ido), 



these new operators satisfy 



^a(^) = (^a(r^))^ d^ = {d^Y; 

{^'„(n), ^/3(n )} = Sa,f3Sn,n', {da, C^/?} = (2) 



Then the Hamiltonian becomes 



3 3 

H = itJ2Yl ^»{n + l)^„(n) + iV ^a(0)ci„ + f/dida^gc/o. (3) 

n a=0 a=0 



When the 0(4) symmetry breaks down to 0(3) symmetry in the hybridization term, the 
model Hamiltonian becomes 

3 

H = itJ2T. ^«(^ + + ^^0^0(0)4 

3 

+lV J2 ^a(0)c?a + Udid2d3do, (4) 
a=l 

where Vq ^ V. In the large U limit (t/ > 0), a Schrieffer- Wolff canonical transformation can 
be applied to generate a s-d type of model, the so-called 'compactified two-channel Kondo 
model', 

^ = E E ^-(^ + + ["^1^(0) + J^m] ■ Sd, (5) 

n a=0 

with Ji = 2V{y + Vq)/U and J2 = 2V(y — Vo)/U, where the local impurity spin couples 
to both the conduction electron spin a(0) and the conduction electron 'isospin' density 
operator r(0) at the impurity site. When Vq = the two exchange couplings are identical, 
i.e., Ji = J2, and it had been conjectured that this form of the model has the same low- 
energy excitations as the two-channel Kondo model ||^. 

To distinguish the Anderson model of Eq.(4) from other Anderson impurity models we 
will refer to it as the 'compactified' Anderson impurity model. Here we concentrate on the 
Vo = case: 

H = Hq + Hj, 

3 3 

^0 = E E ^"(^ + + E "^Mda, 

n a=0 a=l 

Hi = Udid2dzdo. (6) 
This will be used as our model Hamiltonian. The Hamiltonian for Vq 7^ is considered in a 



separate paper 12 



Fourier transforms for the conduction electron operators can be introduced as usual 

1 



-=Y.m^{k)e'^--^ a = 0,1, 2, 3, (7) 
vA u 
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here N is the total number of the sites and the lattice spacing is a. The anticommutation 
relation for the conduction electrons is converted into 



{^a{k),^p{-k')} = 6a,fs6k,k'. (8) 

Substituting these expressions in our model Hamiltonian, up to a constant we get 

^ iV 3 

^ = E E ek^a{-k)^a{k) + E E ^a{k)da + UdMdo, (9) 

fe>Oa=0 ViV 1^ 

where €k — 2tsm{ka) is the dispersion relation of the conduction electrons. 

In the present model, the new 0(3) symmetry in the hybridization is the key feature. 
Since the scalar field {a — 0) decouples from the local impurity, its single-particle Green 
functions, defined in terms of retarded double-time correlation function, is easily found to 
be a free propagator: 

Gl,,{uJn) = - f^dT{TrMk,r)M-k',0))e'^"' = (10) 
Jo lUn - ek 

where Un = (2n -|- 1)tt/13, j3 is the inverse of the temperature. Meanwhile, the vector field 
^a(/c) {a — 1,2,3) hybridizes with the corresponding impurity vector field da, and the 
scattering of the conduction electrons from the local impurity is given by the following 
relation: 

where G^cd^n) is the Fourier transform of the impurity vector propagator defined by 
— {TTda{T)da{T'))H- Then the conduction electron t-matrix is thus expressed as 

tk,-k'{^n) — -j;^Gyec{^n)- (12) 

It is easily seen that this t-matrix only depends on the Green function for the impurity vector 
field Gvec ("?■)■ The implication of this relation is that the conduction electron resistivity will 
be determined by the Gvec('?") only. Moreover, there is also a general relation between the 
impurity vector propagator and the cross correlation function. 



GmK) = {{^a{k)K)) = J^^Z2£(^^ a = 1,2,3. (13) 

This relation becomes useful when we calculate various two-particle correlation functions for 
the conduction electrons, because we can relate the conduction electron correlation functions 
at the impurity site to the impurity correlation functions. 

In addition, we note that under the simple transformation S — V^do, the model Hamil- 
tonian transforms according to 

SH{U)S-^ = H{-U). (14) 

Since the sign change of the on-site Hubbard interaction corresponds to the exchange of 
the impurity charge and spin degrees of freedom, or the exchange of the impurity spin and 
'isospin' density operators, the particle-hole symmetry is kept in the present model, and the 
the chemical potential is pinned at the Fermi level. 



III. THE UNPERTURBED HAMILTONIAN 

Before considering the effects of interactions, it is very constructive to examine the 
unperturbed model Hamiltonian. 

3 3 

Ho^itJ2tl + l)*a(n) + E *a(0)c?a. (15) 

n a=0 a=l 

This is an exactly soluble Hamiltonian in which only the vector field of the conduction 
electrons hybridizes with the impurity vector field, while the scalar field is free. The impurity 
Green functions are defined as follows: 

G'o(t) = -{Trdoir)dom, Ga{r) = -(7;da(T)da(0)). (16) 

It is straight forward to calculate their Fourier counterparts, 

Go{uJn) = T^, Ga{uJn) = —^T , (17) 

where A = TipV"^ is the hybridization width, p = {hvf)~^ is the conduction electron density 
of states, and cun — (2n -|- l)7r//9 is the fermionic Matsubara frequency. Here we find that 



the impurity scalar propagator is a fermionic zero mode with Go{t) = — sgnr/2. Moreover, 
both impurity scalar and vector field propagators are odd in their arguments. 
The local impurity spectral function can be evaluated as 

^'^M = 1:^^X7^ + (18) 

which reveals the basic physics of the unperturbed Hamiltonian that half of the impurity 
degree of freedom is free, while the rest of the degrees of freedom couple to the conduction 
electrons. A similar impurity spectral function was also found in the single-impurity two- 
channel Kondo model IQ. The change of free energy due to the hybridization with the 
impurity can be calculated as 

where /(cij) is Fermi distribution function, and the impurity residual entropy is found to 
be In a/2, revealing that the ground state of Hq is a two-fold degenerate state, which is a 
consequence of the spin-isospin symmetry with [Hq, 5] = for S* = v^rfo- 

In addition, we can calculate the spin and charge density-density correlation functions 
to find out what kind of low-energy behavior is described by Hq. In the Majorana represen- 
tation, the impurity spin and charge density operators are expressed as 

1 i 
Sd = 2^d\d^ - d\di) = --{did2 - 44), 

1 i 
= ^{d\d^ + d\di - 1) = --(44 + 44)- 

The spin and charge density- density correlation functions are equal to 

{T^S^,{r)Sm) = (T.n,(r)n,(0)) = ^ [gI{t) + G„(r)Go(r)] . (20) 
Their corresponding Fourier transform is 

X'^ai^n) = H [Ga{uJn')Ga{uJn - UJn') + G a{uJn')Go{uJn - UJ^')] , (21) 

where uJn = 2mT / (3 is the boson-type Matsubara frequency. The first term in the brackets 
corresponds to the normal FL-like density-density correlation forms, but the second term is 



singular. As far as the singularity is concerned, the imaginary part of the spectral functions 
is obtained by completing the summation over frequency: 

1 A f uj\ cj<r 

ImxTnuJ.T) = -- , , , tanh — ~ <^ ^ 

const., A a; ;» T 

This kind of behavior was assumed in the marginal FL phenomenology to describe the 
normal state of the high temperature superconductors, where the standard FL theory does 
not appear to explain the experimental observations. The conduction electron charge and 
spin density- density correlation functions at the impurity site can also be calculated, but 
their spectral functions do not display any kind of singular behavior. 

It is clear that the Hamiltonian with [7 = has singular behavior due to the degeneracy of 
the ground state and does not describe a non-interacting Fermi liquid. In the perturbation 
theory we develop in the next section, we make an expansion about this unusual weak- 
coupling limit. 

IV. PERTURBATION FORMALISM OF THE MODEL HAMILTONIAN 

A. Partition function and free energy 

Now we consider the full Hamiltonian for the compactified Anderson single-impurity 
model, 

H = Hq + Hi, 

3 3 

^0 = E E ^"(^ + + E *«(o)d„, 

Hi = -Udodid2d3. (22) 
The partition function for this Hamiltonian can be expanded as follows: 



Z = ZoJ2i-ir dTn rfr„_i.... / dn{Hi{T^)Hi{T„,^^)...Hi{n)) 

n=0 ^0 ^0 ^0 

= ZoE^"/ I dTn-1.... (iriF„ (Xn, r„_i , . . . , Ti) , (23) 

n=0 ^0 Jo Jo 
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where Zq denotes the partition function for the unperturbed Hamiltonian, and (...) means 
the thermodynamic average is carried out over the unperturbed part Hq. Noting that 
in Ho the four Majorana fermion components of the local impurity decouple completely, 
FniTn, Tn-i, Ti) cau bc factorizcd as 



((io(r„)(io(r„_i)...(io(ri)) J| ((iQ(r„)(io(r„_i)...do(ri)). 



(24) 



a=l 



Here both the impurity scalar and vector field single-particle Green function Go{t) and 
Gair) {a = 1,2, 3) are free propagators. When the Wick theorem is implemented, it can be 
verified order by order, that each Majorana fermion expectation average can be represented 
by a Pfaffian determinant [0|. For the expectation values of the impurity vector operators, 
the Pfaffian determinant is defined by the square root of an antisymmetric determinant 
composed of the impurity vector propagator G'a(r), 



0, Gairi-T2), Gain-Ts), , Gain - Tn 

Ga{T2-Ti), 0, Ga{T2-r3), , Ga{T2 - Tn 



\D 



Ga{rn - Ti), Ga{Tn - T2), Ga{Tn - T^), , 

As shown in Section II, the impurity vector field propagators are odd in their arguments, 
G a{jr—Ts)+G aiTs—Tj) = and Ga^Tr = Ts) = 0. The Pfaffian determinant \Z}„(ri, T2, .., Tn)\ 
is given by 



Ga{ri-T2), Gain-Ts), Gain - T^) 

Ga{r2-rs), Ga{T2-T„) 

Ga\Tn--l Tfi 



^Ga{Ti-Ta)Ga{Tb-Tc) Gain-Tm), 



where the subscripts l,a,b,c, ,l,m oi each term under the summation are a permutation 

of the first n integers, each Green function Gai^r — t^) has s > r, all different terms of 
this type are included, and the total number of terms is (n — l)(n — 3) • ■ ■ 5 ■ 3 ■ 1. The 
sign attached to each term is positive or negative depending on whether the permutation is 
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even or odd. In addition, the basic property of Pfaffian determinant is that all odd-order 
determinants identically vanish, and all the expansion terms we should consider are even- 
order determinants. For instance, the fourth Pfaffian determinant \D4{ri, T2, t^, T4) \ is given 
by 

I Gain, T2), G 

GaiT2,T3), Ga{r2,T4) 
Gai^s, Ti) 

= Ga{Tl,T2)Ga{T3,T4) - Gain, T-i)Ga{T2-, T^) +Ga{ri, Ti)Ga{T2, T^) . 

On the other hand, the impurity scalar field propagator has a special form, Gq{t) = 
— sgnr/2, its corresponding expectation is trivial to calculate because the imaginary time 
sequence has been assumed P > T2n > T2n-i >■■■■> T2 > ti > 0. 

{do{r2n)do{T2n-i)..Mri)) = (^) ■ (25) 

Therefore, the partition function for the compactified Anderson single-impurity model H 
can be expressed as the cube of the Pfaffian determinant with elements corresponding to 
the impurity vector field single-particle Green function at different times: 

Z / Zq = y--j=^ dT2nJ^ dT2n-l.-.. {\^2n , Ts , . . . , Tg^) | } . (26) 



For the ordinary symmetric Anderson model with Vq = V, the power of the Pfaffian deter- 
minant in the partition function is four rather than three, and in fact there is no need to 



introduce Pfaffian determinant in that case [|1^,|T^]. Thus one power of the Pfaffian determi 



nant in the partition function corresponds to each of the Majorana fermions involved in the 
hybridization. Using the linked cluster theorem, we find that the free energy change due to 
the local impurity can be expressed in power series in the parameter U 



/QN ^ / f/ \ ^" 1 rl^ rr2n rT2 3 

Fimp = F>Jp-22\-^j ^ _/o ^^2n cir2„_l . . . . {\^2n (ti , , . . . , Tsn) | }; , (27) 

where the subscript / on the bracket in the above equation indicates that only linked diagrams 
are to be considered. The extra impurity specific heat can be calculated systematically from 
this free energy expression. 
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B. Single-particle Green functions and self energies 



Now we consider the perturbed impurity single-particle correlation functions. For the 
impurity scalar field do, the perturbed Green function is defined as 

(28) 



Gsc(r,T') 

and its Fourier transform is given by 

p Jo Jo 



{Trdoir)doir')exp 


'-Jo^Hiir)dr]) 


{Tt-cxp 


'-jiHi{r)dr 


> 



(29) 



In a similar way to our early derivation of the partition function, the perturbed scalar field 
propagator can be expanded in powers of U : 

~^ p Jo Jo Jo Jo Jo 
{(T,rfo(r)rfo(r')rfo(r2n).-c?o(ri))\i^2n(Ti,r2,...,r2„)|^}^. (30) 

For the special form of Go{t) = — sgnr/2, the expectation value of the impurity scalar field 
propagator can be calculated as 

-{Trdo{T)do{T')do{T2„)....do{Ti)) 

= E(-l)^^iGo(r - r,)Go(r' - r,) - Go(r - r,)Go(r' - r,)] + [If Go{r - r'). (31) 

Then completing the integrals over r and r', we obtain the scalar field Green function in 
terms of the following equation, 



G'sc(t^n) = GoiuJn) + Go(a;„)S^^(a;„)G'o(a;„), 
where an improper self-energy is represented as 



(32) 



dT2n / dT2n-l- 

Jo 



dTi 



^(-l)*+^[sina;„(r, - r,)] {\D2„(ri,r2, ...,r2„)|}' 



(33) 
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For the impurity vector field da with {a — 1,2,3), the definition of its Green function 
Gyec{^,T') is similar to the scalar field, and its Fourier counterpart Gvec('^n) can also be 
expanded in a power series in U : 

oo / \ 2n ^ , rP n ■ , n 

G^^M = GaK) + E i^-^j ^ I dr dr' dr^n-... dne'^-^^-^ ) 

{{Trda{T)daiT')daiT2n)....dain))\D2n{ri,T2, ...,T2„)|'}^ . (34) 

The remaining expectation average over da can be calculated as 

-{Trda{T)da{T')da{T2n)---da{ri)) 

= E(-l)'^'[G^"(r - n)Ga{T' - Tj) - Ga{T - r,)G„(T' - n)]\Dg\ + Ga{T - r')\i^2n|, (35) 

i<j 

where \£^2nl ^^e so-called cof actor of Pfaffian determinant \D2nl) which is a Pfaffian 
determinant of the same type but with and tj removed from the time sequence. Then we 
can perform the integrals over r and r' and get the impurity vector field Green function, 

G'vec('^n) = GaiuJn) + GaiuJn)Keci'^n)Ga{uJn), (36) 

from which the corresponding improper self-energy is extracted as 

KeM=^^j[^-^j idr^nl dr^n...... I dn 



E(-l)^+'[smcUn(r.-r,)]\D^;(n,T2,...,T2„)|{\D2„(Ti,T2,...,r2„)|}M . (37) 
.i<i )l 

The self-energies usually used are the proper self energies, which are related to the 
improper ones by 

Ssc(<^n) = Ssc('^n)[l +Go(u;n)Ss^(<^n)]"^ 

S.ec(t^n) = KeM[l + G a{uJ n)Kec{^ nT^ ■ (38) 

As far as the leading-order perturbative contributions are concerned, it is not necessary to 
distinguish the improper self energies from the proper ones. 
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C. Two-particle Green functions and vertex functions 

The method of evaluating the single-particle correlation functions can be applied further 
to the impurity two-particle Green functions. Then the two-particle vertex functions can be 
calculated to give a more complete picture of the physical behavior of this model. Generally, 
three different two-particle correlation functions can be defined. The first one is given by 

^ ... f dTdT'dT''dT''\T,do{T)d^{T')d2{T'')d^{T'''))e^ (39) 

Jo Jo 

Using the same strategies we used earlier in calculating the impurity single-particle Green 
functions, we find that 

= -x/2GoHG„(a;0GaMGaK') E [^J i '^^'-+^ i dr^n-... dn 
The corresponding improper vertex function can be defined as 

°° \/2 f U \^"+^ r/S rT2n+l fT2 

r[,,,,,3(^, u\ u\ J") = - E y (^^ j I dr,^^, I dr,^.... I dn 

I »J ) I 

When all the frequencies are set to zero, this expression simplifies to 

r(),i,2,3(o, 0,0,0) 

= E V (^j I {i:(-l)-\«k«(n....,r.„«)l[ . (42) 

This vertex function will play an important role when we apply the multiplicative renormal- 
ization group to our perturbation expansion in section VII. 

The above treatment can be applied to the following two two-particle Green functions 
straight forwardly: 

{Trd^{T)do{T')da{T")do^{T"')) and {Trda{T)d^{T')dp{T")dp{T"')), 
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here a, P — 1, 2, 3, and (3. Their resulting equations of motion turn out to be 

= P^G,{u)G^{u")5^,,_J^,.,_^. + /3G'o(^)G'o(u;')G'aMG'.(a;'")r[,,o,a,a(c^, ^"') ; 

where the improper vertex functions are given by 

gi(u;riWr,) _ gi(u;r,Wr,)j Jgi(u;"v Wr,,) _ gi(c."r,, |\D2„| }^ . (45) 

These two improper vertex functions are antisymmetric in their arguments and 

r[),o,a,a(0, 0, u", J") = n,o,a,a(^, 0, 0) = r:,,„,^,^(0, 0, o;", u'") = r:,,,,^,^(^, u;', 0, 0) = 0. 

Since the last two vertex functions vanish in the zero frequency limit, we expect only 
To,!, 2,3 to be important in determining the two-particle interactions in the low-energy regime. 
It is straightforward to relate the improper and proper vertex function, for example, 

10,1,2,3(0;, u;', a;", a;"') = ^^T',^^^^^^{u,u\u\uj"'), (46) 

where Zsc{oJn) — Gscii^n)/Go{cUn) and 2;vec('^n) — Gvec{u!n) /Ga{uJn) denote the wave function 
renormalization factors for the impurity scalar and vector components, respectively. 

D. Dyson equations of the impurity self-energies 

The impurity scalar and vector self energies Egc and Eyec have already been derived in 
terms of the Pfaffian determinants and their co-factors, which are functions of the unper- 
turbed propagators. However, when it is necessary to sum an infinite series, instead of just 
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calculating the first few diagrams, it is usually more convenient to express the impurity self 
energies in terms of the vertex function Fq, 1,2,3- On the other hand, these two impurity self 
energies are also dependent of the vertex function ro,i,2,3- The procedures are illustrated by 
considering Sgc as an example. 

Since Esc is even function of U, there are no first-order perturbation corrections. Also 
there are no higher-order diagrams in which the self energy is connected to the basic prop- 
agator Gsc by a single first-order vertex function Fo]] 2,3(^^1 ^"j ^^'") = ~^^a;+w' w+a;"',o- 
The second-order perturbative corrections to the self energy Sgc will be derived in Section 
VII, here we just use the result, 

E(2)(cj„) = -— ^ G«(cui)G«(cj2)G^K - ^1 - ^2). (47) 

It is not hard to see that all the higher order diagrams for Sgc can be obtained from this 
second-order diagram by insertion of self energy parts into internal unperturbed propagators 
Ga and replacement of the square on the right by the proper vertex function ro,i,2,3- This 
diagram is displayed in Fig. la. Then the Dyson equation of the impurity scalar self energy 
is derived in terms of the impurity vertex function, 

= G'vec(t^l)Gvec(t^2)G'vec(t^n - UJi - U;2)Fo,l,2,3(t^2, UJn - UJi - UJ2; U)l, ~UJn). (48) 

In a similar way, we can build up the Dyson equation for the impurity vector self energy 
Svec from its second order perturbative corrections, which will also be derived in section VII, 

SE(^n) = -S E Go.iuJ,)G^ic02)GoicOn ' UJ, - U2) . (49) 

UJl,UJ2 

The corresponding Dyson equation is given by 

= G^cc{^^l)G^cc{uJ2)Gsc{uJn - LJi - t^2)Fo,l,2,3(^2, - C^i - UJ2; UJi, -Un), (50) 

the corresponding diagram is described by Fig. lb. Eq.(^8D and Eq.(^) are two basic equa- 
tions of the model Hamiltonian and also exhibit the relations between the impurity self 
energies and the two-particle vertex function which is finite in the zero frequency limit. 
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V. DYNAMICAL SUSCEPTIBILITIES 
A. Impurity spin and charge dynamical susceptibilities 

As noted earlier the impurity charge and spin density operators are defined by 

Sd ^ --{did2 - dods), --{did2 + dods). (51) 

Using these expression the impurity spin and charge density-density correlation functions in 
the Fourier space can be expressed in the form, 

2 



XT'^M = ^{XevenM " XoddM}, (52) 



with 



XevenM = \ f dr ['^ dr' {{Trd,{T)d,{T')d2{T)d2{T')) + {Trdo{T)do{T')ds{T)ds{T'))} , 

p Jo Jo 

Xodd(^n) -IfdT r dT'{T,do{T)d,{T')d2{r')ds{T)). (53) 

p Jo Jo 

Expressions for these correlation functions in terms of the Pfaffian determinants are 
rather complicated. We can express them however in terms of the impurity improper self- 
energies and vertex functions, which is a more convenient form for the calculation of the 
low-order perturbation corrections. For the even part, the form is given by 

Xeven(^n) 

112 

^ ■^Z]<^«('^i)Ga(a;„ -a;i) + -J^GaMGoiuJn - ^^i) - -^J2'^lM'^<^i'^ri - '^l)KecM 

r-' U)i r-' tJi 

E GaMGaMGaiuJn - UJi)GaiuJn + U;2)r^,a,/3,/3(^l> ^^n " ^1, "^n " ^^2) 

P UJl,U2 

E Go{uJi)GQ{uJ2)Ga{uJn - i^l)Ga{uJ„ + UJ2)rQ^0^l3^/3{uJi,UJ2;UJn - UJi, -UJn - (54) 

P U)\,U}2 

while the odd part is given by 
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Xodd(<^n) 

2 

^ Yl Go(ct;i)Ga(cU2)Ga(ti^n -'^l)GQ(c^n + ^2)ro_i^2,3(^l'^2;a;n -CUi,-CUn -a;2), (55) 

where Go and Ga correspond to the unperturbed impurity scalar and vector propagators, 
respectively. 

If the bosonic Matsubara frequency cUn is taken to zero, the static charge and spin 
susceptibilities are obtained. In fact, there is another way to derive the impurity static 
susceptibilities. When an additional term SH is added to the Hamiltonian corresponding 
to a coupling with an external field, such as a coupling to a uniform magnetic field 5Hj — 
—i{gfiBh/2){did2 — d^ds) or a chemical potential term 5H — —i{ij,/2){did2 + dod^), the 
partition function can be written in the form 



r Hi{T)dT exp - [%Hi{T')dT' 
Jo Jo 



(56) 



Then one can expand this partition function to second-order in SH in a power series in U, 
and derive expressions for the static spin and charge susceptibilities: 

imp f 9f^B \ J ~ I ~ \ imp / ~ ~ 1 

Xa I 2 J iXeven ~r Xodd / j Xp ^ \Xeven Xodd/ j 

where Xeven and Xodd a-re expressed in terms of the Pfaffian determinant and its co-factors, 
respectively. Since we do not use these forms to calculate the static susceptibilities, we shall 

not write them here. 

It should be pointed out that the above expressions for Xeven (<^n) and Xodd('^n) can also 
be converted into alternative forms in terms of the perturbed propagators Ggc and Gvec, and 
the proper vertex functions Ta,a,/3,i3 and ro,o,a,a, which are general expressions for the model 
Hamiltonian. 

Xeven ((^n) = X) '^"('^l)<^a('^n " ^^l) + ^ X] '^"('^l)<^o('^n " ^^l) 

2 11 

-^^G^ec{^l)Ga{uJn - ^l) - -^^Gyec{^i)Go{Un - UJi) - -5 ^ G'sc(t^l)G'a(a;„ - UJi) 

2 11 

+ -5 X/ ^vec('^l)G^vec('^n ~ ^l) + "5 X/ ^vec('-^l)Gsc(i-i-'n — ^l) + ^ Gs,c{^l)G^ec{^n — ^l) 

P U>\ P P U)l 
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~1q2 X/ ^sc(^l)G'sc('^2)G'vec('^r!, — ^l)G'vec('-^n + ^2)'r0fl,f3,l3{^l, ^2] ~ ^l, ~'^n ~ 1^2); 
ljJl,u>2 

Xodd('^n) 
2 

= XI Gsc(cc;i)G'vec(t^2)Gvec(c^n - UJi)G^ec{uJn + tU2)ro,l,2,3 (t^l , ^^2; tU„ - LOi, -UJn - UJ2) ■ (57) 

These two expressions can be described diagrammatically in Fig. 2, where each term corre- 
sponds to one Feynman diagram. 

By the way, the free energy (pTf) can also be converted into a function of the impurity 
perturbed propagators and vertex function Xodd(i^n) when we use the relation 

\D2n(ri,r2,...,r2n)| = -;^(-l)*G'„(r2„-ri)\D^„_i|, 

i 

and then integrate out the imaginary time Tm The result is 

1 ^ 

{F - Fo)imp = ^Y.J Xodd{uJn)dU 

vcc 

ro,l,2,3(t^l, ^^2; t^n - tUl, -t^n " ^^2) ■ (58) 



B. Local conduction electron dynamical susceptibilities and singlet superconducting 

pairing correlation function 

The conduction electron charge and spin density operators are defined in terms of the 
corresponding Major ana fermions, 



-- [*i(n)*2(n) - '^o{n)'^3in)] ; 

i 
2 



(59) 



S'An) = 1 [G\in)C^in) - C|(n)Q(n); 

ndn) = ^ [c\{n)C^in) + Cj(n)Q(n) - l] = [^i(n)vi>2(n) + <iIoin)^3in)] 
and their density- density correlation functions at the impurity site are given by 

{TrS',{0,T)S',{0,T')) = (T.ne(0,r)n,(0,r')) 

= ^ {(T,^i(0, r)^i(0, r')^2(0, r)^2(0, r')) + (T.^o(0, r)^o(0, r')) (T.^3(0, r)^3(0, r'))} , (60) 
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where the conduction electron scalar operators decouple from the vector operators because 
they are free Majorana fermions in the model. Due to the hybridization between the impurity 
vector field and the conduction electron vector field, the spectral function of these two 
density-density correlations can be derived from the perturbation expansion series of U in 
the same way we did for the impurity density-density correlation functions. The final result 
is expressed as 

E Gki,d{t^'l)Gk2,d{^2)Gk3,d{^n — ^l)Gk4„d{^n + '^2)ri,l,2,2('^l> "^2; — <^1; —^n — ^2) 
P UJ\ ,ijJ2 

-■^ E G'fci,d(a;i)Gfc2,rf(^i)G'fc3,fe4K - ^i)^vec(^i) 

-^T.Gks4MGk,,dM [G°,„fe,(a;„ - a;i) + G^^^^iuj^ - uj{)] ■ (61) 

Meanwhile, we can consider the singlet superconducting pairing correlation function at 
the impurity site which is defined by 

X'ZM = ^ /^""cir j^''rfr'(T.Q(0,r)Q(0,r)C|(0,r')C|(0,r'))e^-"(— (62) 

where a;„ = 2m:/ (5. Using the Majorana fermion representation, we obtain the following 
expression 

(T,Q(0,r)Q(0,r)C}(0,r')C{(0,r')) 

= ^ {(7;*i(0, t)*i(0, t')*3(0, t)*3(0, t')) + (7;*2(0, t)^'2(0, t')*3(0, t)*3(0, t')) 

+ (T,^o(0, r)vl/o(0, r')) (T,vl/,(0, r)vl/i(0, r')) + (T,^o(0, r)v]/o(0, t')){T,^^{Q, r)^2(0, r'))} , (63) 

and, a relation with is found 

X:up(^n) = 2x^:?(c.n). (64) 

For the present Anderson-type model, we have several general relations between the 
conduction electron propagators and the impurity propagators. These relations for the 
Majorana Green functions are 
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Gk^A^V = -Trf- ■ (65) 



Substituting these relations into the expression (|61|), we finally get the result 

■(vrp)^ ^ G'Q(u;i)G'a(c^2)G„(t^„ - u;i)Ga(a;n + u;2) 



r'l, 1,2,2(^1^2; - cui, -t^n - cj2)sgncjisgncj2sgn(cjn - u;i)sgn(t^n + ^^2), (66) 

where the density- density spectra in the case of f/ = are 

1 W"^ 
= --5(7rp)^^sgnajisgn(aj„ - Ui) + z— — (yrp)^ ^ Gc.(wi)sgn(w„ - ^1) 

LUX P 1^1 

+ ^ (^P) ' E G'a (u;i ) G - 071 ) , (67) 

^P OJl 

which is a regular contribution. When the Matsubara frequency is taken to zero, the cor- 
responding static susceptibilities are also obtained. With the help of the relation (p3), the 



conduction electron singlet superconducting pairing correlation function at the impurity site 
can also be obtained. 



VI. MAIN PERTURBATIONAL RESULTS OF THE MODEL 

A. Leading specific heat correction 

In the section IV, an expression for the impurity free energy was derived in a power series 
in U . Since F^^^^ is the corresponding free energy in the case of f/ = and is just a regular 
contribution, the first singular contribution to the free energy comes in second order in f/, 
and is 
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(2) U'^ 

^imp = -Jn dT2 dTiGHn - T2)sgn{T2 - Ti) 



IP [0 

^ X) Ga{uJl)Ga{uJ2)GQ{uj^)Ga{uJl+UJ2+UJ2). (68) 



The singular contribution from this term is 



which gives a singular correction to the impurity specific heat, 



B. Self energies and electrical resistivity 

Prom the general expressions for the self-energies, it can be seen that the lowest order 
self energy corrections are of second order. For the impurity scalar field, the perturbed self 

energy is 

T.fX^n) = -2^^ f dT2 r dnGlin - ra) sin[a;„(ri - T2)] 
p Jo Jo 

f/2 rf' 



= -— / Gl{T)e''^^^dT. (69) 
2 J — [3 

In the imaginary time space, this self energy is 

^'S{t) = -U'GI{tI (70) 

which can be represented by a Feynman diagram shown in Fig. 3a, and its Fourier transform 
is given by 

E(2)(^„) = --55- E G^{uj^)Ga{uj2)G^{uJn - ^1 - UJ2). (71) 



The imaginary part of its retarded Fourier transform after integration is given by 

(72) 



ImE2'(..r)»4(^)^ 



a) + (x) 
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Here the analytical continuation has been performed, and according to the Kramers-Kronig 
relation the real part should be 



ReE(^)(a;,T)«.-(-^)'a;. 



(73) 



Moreover, the renormalization factor of the wave function or quasiparticle weight (2;sc('^n) 
taken at a;„ = 0) to second order is given by 



.(2) 



a)=0 



(74) 



The scalar self energy exhibits the normal local FL behavior, and such a behavior ensures 
the fermionic zero mode in the unperturbed Hamiltonian is preserved in the presence of 
perturbations. 

On the other hand, the vector perturbed self energy can be written down from the general 
expression. 



^vJci^n) = -i-T / d^2 / dTiGKn - T2) sin[u;„(Ti - T2)] 
p Jo Jo 



(75) 



The sign function appears when the double integrals are reduced to a single integral. In the 
same way, the self energy in imaginary time space becomes 



Stl(r) = -U'Gl{r)Go{r), 



(76) 



which corresponds to the Feynman diagram in Fig. 3b, and the Fourier transform is given by 



^lec(^n) = GaMGaMGoiuJn - UJi - UJ2). 



(77) 



Since this self energy involves the fermionic zero mode, the calculations of its spectral func- 
tion becomes subtle. The retarded imaginary part is given by 



A»|a;|»T. 



(78) 
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The corresponding real part is also obtained as 



(79) 



where x — max(|a;|, T). Such a self energy is greatly different from the form given by the 
ordinary FL theory, in particular, the renormalization factor of the wave function or the 
quasiparticle weight logarithmically approaches to zero as T — > 0, 



vec 



do; 



u>=0 



H -1 



(80) 



which implies that the impurity vector self-energy displays the marginal FL behavior. 

There is a singular temperature-dependent contribution to ImE^g),(0, T), and the conduc- 
tion electron t-matrix determined by the Green function of the impurity vector propagator 
has been derived as follows 



N 



Gyec{l^, T). 



(81) 



Thus, the imaginary part of this t-matrix can be expressed in terms of the retarded self 
energy for the impurity vector field: 



Imi(a;, T) 
1 



+ 



ImS^ 



(82) 



T:pN {u - ReSvec)^ + (A - ImE^^^f npN {cu - ReEvec)^ + (A - ImEvec)^ ' 
where the first term describes the elastic scatterings, while the second term describes the 
inelastic scatterings. 

Assuming the conduction electrons incoherently scatter from the dilute magnetic impu- 
rities, the linear response theory allows the electrical conductivity to be expressed in terms 
of the following Kubo formula 



(83) 



here f{x) is the Fermi distribution function, Vf is the Fermi velocity of the conduction 
electrons with charge e and density of states p, and t{u!,T) is the electron relaxation time. 
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which is related to the t-matrix 



T 



-1 



2Ni^pImt{u;,T) 



where A^imp is the total number of the impurities, and riimp = Ni^ap/N is the impurity 
concentration and is supposed to be much less than unity. On substituting the second order 
self energy for the impurity vector field into equations (82) and (83), we obtain the electrical 
resistivity in the low temperature regime. 



Here we have just taken into account the contributions of order (:;^)^- Therefore, up to 
the second order of perturbation series, the electrical resistivity has a linear temperature 
dependence, which is a direct consequence of the anomalous behavior of the imaginary part 
of the impurity vector fermion self energy. Such a resistivity makes the weak coupling fixed 
point of the present model differ from the strong coupling fixed point of the two-channel 
Kondo model. 



As we said earlier, in the zero-frequency limit there is only one non-zero two-particle 
vertex function Tq^^ 2,3; the other two vertex functions vanish. Its general expression in 
terms of U has been given in section IV. In the first order perturbation, it is trivial to obtain 
as 




(84) 



C. Vertex functions in the zero- frequency limit 




(85) 



the corresponding diagram is described in Fig.4a. To order U^, it is given by 



rS2,3(o, 0,0,0) 
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Using the diagrammatic method, we can clarify that four basic Feynman diagrams have 
contributions to this zero-frequency vertex function, and then re-group these terms to get 
the following expression. 



-GI(ti - T2)Ga{n - T3)sgn(ri - T2) - ^^(ti - T2)Ga{n - n)sgn{T2 - Ts) 



rg2,3(0, 0, 0, 0) = ^ j'^ dTs dT2 dniZGlin - T2)G^{n - r3)sgn(ri - T3) 



-Gl{n - T2)G'„(ri - r3)G',(T2 - r3)sgn(Ti - T2)}, (87) 

where each term is described by the respective Feynman diagram (b), (c), (d), (e) in Fig.4, 
and it is easily found that only the first term (Fig.4b) gives the leading singular contributions, 
while the remaining three terms can be neglected in the hmit T — > 0. Thus, the leading 
third order perturbation correction to the vertex function is 



-^^G{ljJn)Go{uJn) 



rg2,3(o,o,o,o) = -3t/^ 

«^-3t/(^)'ln(^), for T«A, (88) 

and a logarithmic singularity appears. Up to third order in t/, the proper vertex function is 
calculated and exhibits the same singularity as well, 

ro,i,2,3(0, 0, 0, 0) = ^-^^-e',r[),i,2,3(0, 0, 0, 0) « -[/ |l + (^)' + 6 (^)' in (^) I , (89) 

which implies that the higher order perturbative expansion terms will become as important 
as low-order contributions in determining the low-temperature behavior of the model. At 
least the leading order logarithmic terms have to be summed to obtain the correct low- 
temperature behavior. The summation of all the leading logarithmic contributions for the 
vertex function is the so-called parquet approximation. In section VII, we will use another 
equivalent approach: the multiplicative renormalization-group method, to find the final 
results. 
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D. Impurity static susceptibilities 



We have already obtained general expressions for the impurity dynamical susceptibilities 
in terms of the impurity scalar and vector self energies and their two-particle vertex functions. 
Since these associated self energies and vertex functions have been obtained in the lowest 
order perturbation theory, we can use them to derive the even and odd static susceptibilities. 
In the zeroth order, there are two contributions. 



y{0) 

A even 



1 



+ 



1 



(90) 



where the first term is a regular FL like contribution, while the second term shown in Fig.Sa 
is singular, and can be expressed as 



1 



(tta) 



1 



A \ 



(91) 



where ip{x) is the digamma function, and when T <^ A, (| + ~ '4^ (i) ~ (^)' 
and —^Y^ui^G'^ai'^n) ~ {j^- Thus, the singular contribution gives rise to a logarithmic 
divergence in the impurity static susceptibilities. In the first order in [/, there is also one 
singular term corresponding to the diagram in Fig. 5b. From the first-order vertex correction 
for ^'o,i,2,z{'^i^ ^2; -t^i, -^2), we get 



v^^^ -2U 

Xodd — ^'J 



(92) 



While to second order in [/, the perturbative corrections to the impurity static susceptibilities 
become complicated. 



X 



1 

7 



2 r 



1 



-^G{u;n)Go{u;n) 

P <^n 

-^Y^GMGoM 

P U!n 



+ less singular terms. 



(93) 



where we have used the second-order results for the impurity self-energies E'4^^(a;„) and 

^veci^n) and vertex functions r'/^^''2_2(^i5 ^2; — '^i, —^2) and Fg^Q^g 3(tc;i, tc;2; — —uj2)- Note 
that apart from the logarithmic contributions represented by the diagram in Fig. 5c, there 
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appears another singular term: a squared logarithmic term arising from diagram (d) in 
Fig. 5. Therefore, up to the second order in the perturbation series, the singular parts of the 
impurity spin and charge static susceptibilities are given by 



\ 2 J 
1 



A 



A 



(^) (^)^"\t 



VttA 



(tta) 



A 



Vtta; \t 



1 



A 



ri-vr-)in% 



(94) 



E. Conduction electron local static susceptibilities and singlet superconducting 

pairing susceptibility 

For the local static susceptibilities of the conduction electrons, it is straight forward to 
obtain results for the local static susceptibihties of the conduction electrons by taking the 
zero-frequency limit for the dynamical susceptibilities. As we have stated before, the un- 
perturbed charge and spin static susceptibilities of the conduction electrons at the impurity 
site are regular. Then, after subtracting the contribution of the free conduction electrons 
(y = 0), it is expressed as 



1 



(95) 



Unlike the impurity static susceptibilities, the first perturbative correction comes from sec- 
ond order in U. Since the vertex correction Tii2.2{^i:^'2',^^n — ^ii—^n — ^2) involves a 
singular contribution analogous to the diagram Fig. 5c, the singular part of the second-order 
correction is 



con.(2) ^ -^{T:pVY 



1 



1 



-^Go,{Un)Go{Un) 



(96) 



When T ^ A, up to the second order in [/, the static susceptibilities are given by 



(97) 



where there appears a logarithmically temperature dependent contribution. The possible 
reason is that the impurity vector field transfers its low-energy singularity into the local 
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conduction electrons through the hybridization, so there is also a logarithmic temperature 
dependence of the impurity static susceptibilities to second order in U. 

Another important result is the second order correction to the singlet superconducting 
pairing correlation function at the impurity site, its static susceptibility can also be obtained 
through the relation xtnp — ^Xp^J?; 




which indicates that the conduction electrons form a singlet pairing resonant state around 
the local impurity site, different from a single-particle resonant state in the single- channel 
Kondo problem. It should be pointed out that this pairing resonance is pinned at the Fermi 
level (chemical potential). 



VII. MULTIPLICATIVE RENORMALIZATION-GROUP ANALYSIS 

Since the perturbation theory on the present single-impurity model gives logarithmic 
contributions to the impurity vertex function and vector field propagator, the summation 
of the leading order logarithmic terms should lead to some form of scaling behavior. One 
way to sum these terms is by using the usual parquet approximation. Here we use another 
equivalent method to extend our second order perturbational results for the impurity self 
energies and vertex function: the multiplicative renormalization-group (RG) 

The usual multiplicative RG is a simple transformation procedure in which the Green 
functions, vertices and coupling constants are multiplied by real, frequency independent 
factors. The requirement that the two-particle interaction form defined in the model Hamil- 
tonian be satisfied by both the original and the transformed quantities gives a relation 
between them. The classical formulation was used to obtain scaling laws for the X-ray edge 



problem and the single-channel Kondo problem |T7|JT8| . For the logarithmic problem, gen- 
erally, it is important to give a proper definition of the invariant coupling, the temperature 
dependence of which characterizes the behavior of the system. On the physical ground, we 
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take the hybridization width A as a natural scaling parameter. The perturbed impurity 
Green functions and the important two-particle vertex function are written in the form, 



c/f 0,1,2,3 ^ - (^) ro,i,2,3(o, 0, 0, 0), 



(99) 



where Go and Ga are the unperturbed Green functions and tJ is the dimensionless bare 
vertex. For simplicity, in the functions of Zgc, Zy^c, and ro,i,2,3, the frequency variables are 
fixed at the Fermi level and only the temperature variables are retained. If the interaction 
is cut off at the energy A, the Green functions and vertex depend, as a rule, on the relative 
energies A/T. 

Multiplicative RG is usually formulated as the transformation induced by a change of 
the cut-off parameter from A to A', i.e.. 

/A' 



\ T 

^vec I ) U' 



Z2Zvec [—,U], 



ro,l,2,3 f-^, U'j — Z^ ^ro,i,2,3 (^j;, u ] , 



U' 



-1/2 -3/2 j-j 
Zi Z2 Z3U , 



(100) 



where zi, Z2, and Z3 are independent of the temperature variables. Such a transformation 
can keep the four-legged vertex unchanged, which corresponds to the two-particle interaction 
of the model Hamiltonian. Since z^l^ = [1 + U^]~^ does not depends on the cut-off factor 
at all, it can not induce any essential renormalization. For simplicity, we choose Zi = 1. 
Therefore, when the above relations are obeyed, the cut-off dependent U', a self-consistent 
solution to these equations is given by 



vec I ~rj^ 1 



3/2 



= c/ro,i,2,3(^,f/ 



A - 

vcc I T^Ti U 



3/2 



, (101) 



or 
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U' 



^^0,l,2,3(f , U) 



^vec( y ) 



3/2 



3/2 • 



(102) 



ro,l,2,3( 7^, U') 

which is the invariant couphng constant of this model. When the denominator is normahzed 
to unity at T = A', a simple definition of the invariant coupling constant is 



C/inv = C/ro,l,2,3(|^,^ 



u 



1 3/2 



(103) 



Up to the second order in U, the impurity renormalization factors of the wave functions 
and the vertex function have already been obtained as 

1 



^VeC [ rpl '-^ 



To,!, 2,3 ( jT) U 



l + C/21n(f)' 



(104) 



Actually, the non-logarithmic parts are not essential in the present RG approach. Prom the 
above second-order results, the multiplicative factors Z2 and can be derived 

1 1 



Z2 = 



1 + C721n(f)' i + 6C721n(f)' 

and then the dimensionless invariant coupling is given by 



(105) 



(106) 



from which the basic RG scaling equation can be cast into a differential form. 



d\nU 



dlnA 2 



(107) 



The characteristic feature of this scaling equation is that the dimensionless coupling constant 
U = U / (tt A) always increases as the high-energy scale A is reduced. Integrating this 
differential equation gives the relation, 

[7^ 



U' 



(108) 



l + 9C72 1n(f)' 

The scaling relation can be used to extend the second-order perturbational results, and 
effectively sum the leading order logarithmic terms. The scaling trajectory can be traversed 



32 



from the initial parameters A and U, to an effective A' of order T, and an effective coupling 
constant U', which becomes temperature dependent, 

U\T) 



1 + 9t/2 ln(j^ 

When we substitute this new temperature dependent coupling constant into the second- 
order perturbational results, we will get the results which are equivalent to a summation 
of the leading order logarithmic terms. For example, the electrical resistivity in the second 
order perturbation theory was calculated to be 



imp 



(109) 



After substituting the new couphng parameter, it becomes 



P(^) = — ^ 



1 + 



(110) 



l + 9f/2 1n(j) V A 

The In T term is in the denominator due to the summation of the leading order logarithmic 
series. As T — 0, the difficulties arising from the InT term become more severe as there is 
a divergence at finite temperature. 



Tc = Aexp 









if)] 



;iii) 



which is a new weak-coupling low-temperature energy scale. When T > the perturba- 
tional scaling can be extended down to an effective hybridization width A' and effective 
coupling constant U', while T < T^, the electrical resistivity derived from perturbation 
scaling diverges. 



In the conventional RG treatments for the single-impurity Kondo problems |]T9|, the 
conduction electron band-width D is usually chosen as the high energy cut-off factor, and 
the coupling parameters are renormalized as the bandwidth D is decreased. However, for the 
Anderson-type impurity model, the perturbation expansions in the interaction parameter U 
by Yamada and Yosida |T^JT^ and the numerical RG calculations have shown that the 
model provides its own intrinsic high-energy cut-off of the order of the hybridization width 
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A (A > U), inconsistent with using the conduction electron bandwidth as the effective 
cut-off factor . Similarly, the present perturbation results of the compactified Anderson 
impurity model have also proved that the parameter A plays the role of the high-energy cut- 
off factor. In fact, the behavior of the present model is independent of the conduction electron 
bandwidth, which can be taken to the infinite bandwidth limit without lost of generality. It 
is noteworthy to point out that the scaling equation ( |107|) is a for the dimensionless coupling 
parameter U. Due to U = U/{ttA), there is another form of the RG scaling in terms of 



(tta) 



. However, this form of the 



the dimensional coupling parameter U, -jj^ = U 
scaling equation is not useful for discussing the scaling behavior of the interaction because 
the parameter U is not the RG invariant coupling parameter. The RG scaling of the model 
has exphcitly displayed in Eq. (|107|) . 



What have we learned from these scaling arguments? First of all, the hybridization 
width A can be reduced dramatically for the calculation of thermodynamic behavior from 
its initial value down to the thermal energy scale and still be described by a model of the 
same form but with a renormalized coupling constant. Second, there is a subtle point in 
the perturbation treatments. Basically the model has been described by two independent 
coupling parameters A and U, and the perturbation expansion is based on the following 
two-particle vertex function, 

ro,l,2,3(0,0,0,0) = -f/fo,i,2,3(f/), (112) 

where ro,i,2,3 is the corresponding dimensionless vertex function and only depends on the 
dimensionless coupling constant U = The important point is that dimensionless coupling 
parameter \U\ increases as the scaling parameter is reduced. U becomes infinite when the 
scaling parameter approaches the characteristic low-energy scale Tc, and the dimensionless 
vertex function diverges as well. Thus, the perturbation theory begins to break down, and 
there is one weak-coupling low-energy scale depending on the dimensionless coupling 
parameter U only. 
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VIII. DISCUSSIONS AND CONCLUSIONS 



From the previous perturbational calculations and the multiplicative RG analysis, we 
have found a non-FL behavior in the weak-coupling regime ([/ -C 1), and that the invariant 
coupling parameter U increases as the high-energy scale A is reduced. So our perturbation 
results are not valid when the coupling constant becomes larger and larger, or the energy 
scale is decreased to the low-temperature regime (T = T^). 

However, in the large-?/ limit, the Schrieffer- Wolff canonical transformation can also be 
applied to the present model, and a s-d type of model (so-called compactified two-channel 
Kondo model [H) is obtained 

3 

H = ^tY.Y. ^"(^ + l)^a(ri) + J [a(0) + r(0)] ■ 5^, (113) 

n a=Q 

where J = 2V'^/U and pJ = (^)' weak-coupling limit pJ ^ 1 of this s-d model 

corresponds to the large-?/ limit, and the multiplicative RG analysis can be also used to this 
model as the usual treatments for the ordinary Kondo model [|T^|18|. The scaling equation 
has been derived as 

d{.pJ) _ ro/„ T^2 / TN3 



-[2{pJY - {pj)% (114) 



where the conduction electron bandwidth D plays the role of the high-energy cut-off. 
This scaling equation explicitly shows that the dimensionless coupling parameter pj grows 
stronger and stronger under the RG transformations. For the corresponding compactified 
Anderson model in the large U regime, this is equivalent to a decrease in the dimensionless 
coupling parameter U {pJ = 2/ (vr^f/)) as the low-energy scale is reduced. Using the relevant 
cut-off factors in the small and large U regime, we arrive at a schematic flow diagram as 
shown in Fig.6, and we conjecture a stable fixed point of the model Hamiltonian in the 
intermediate coupling regime at U = Uc- This flow diagram is analogous to that of the 
isotropic two-channel Kondo model Obviously, there are significant differences between 
the present compactified Anderson impurity model and the two-channel Kondo model. 
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In conclusion, we have developed a systematic perturbation theory for the compactified 
Anderson impurity model in the Majorana fermion representation, where the unperturbed 
Hamiltonian has a degenerate ground state. We have calculated the leading perturbational 
corrections in the weak coupling limit. We have also derived some general relations for 
the impurity susceptibilities and self energies in terms of the vertex functions. The main 
results of the paper are that a linear temperature dependence of the electrical resistivity is 
obtained from the second order theory, and some non-FL thermodynamic properties have 
been calculated as well. The conduction electron singlet superconducting pairing correlation 
function at the impurity site is found singular in second order in U, indicating the formation 
of the singlet conduction electron pairing resonance at the Fermi level. In the third order in 
U, the vertex function ro,i,2,3(0, 0, 0, 0) has logarithmic corrections, and the summation of 
the leading order logarithmic terms results in a new weak-coupling low-temperature energy 



down. The behavior of the low-energy excitations below Tc is still open for the future 
investigations. 
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below which the perturbational approach begins to break 
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FIGURES 

FIG. 1. The Dyson's equations of the impurity self energies in terms of the impurity vertex 

function ro,i,2,3- (a) is for Sgc, and (b) is for Svec- The solid lines denote the perturbed propagators 
Gvec and the dotted lines are Ggc- 

FIG. 2. The even and odd parts of the impurity spin and charge dynamical susceptibilities in 
terms of the impurity propagators and vertex functions. The thick solid lines denote the perturbed 
propagators Gvec and the thick dotted lines are Ggc, while the thin solid lines correspond to the 
unperturbed propagators Ga and the thin dotted lines to Gq. 

FIG. 3. The second-order perturbative corrections to the impurity self energies, (a) is for the 
scalar field, and (b) for the vector field. 

FIG. 4. The first order and third order corrections to the impurity vertex function 
Tq 1 2 3(0)0)0)0)- (a) is the first order correction, (b), (c), (d), and (e) are the third order contri- 
butions. But only the diagram (b) gives the logarithmically temperature dependent correction. 

FIG. 5. The perturbative corrections to the even and odd part of the impurity spin and charge 
static susceptibilities in the low order perturbations, (a) is the singular diagram of the even part in 
the zero order in U. (b) is the odd part in the first order perturbations, (c) is the In T contribution 
to the even part in the second-order of U. (d) is the diagram giving rise to the In^ T terms of the 
even part in the second-order perturbations. 

FIG. 6. The schematic RG flow diagram in terms of the dimensionless coupling parameter 
U = U/{TrA). The arrows denote the direction of decreasing the relevant high-energy scale. Uc 
corresponds to the stable fixed point of the model. 
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